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Abstract: MHD effects on peristaltic flow of a 
Bingham fluid in a channel with permeable walls is 
studied under long wavelength and low Reynolds 
number assumptions. This model can be applied to 
the blood flow in the sense that erythrocytes region 
and the plasma regions may be described as plug 
flow and non-plug flow regions. The effect of yields 
stress, Magnetic parameter, Darcy number and slip 
parameter on the pumping characteristics is 
discussed through graphs. 
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I. INTRODUCTION 

Physiological fluids in animal and human bodies are in 
general, pumped by the continuous contraction and 
expansion of the ducts. These contractions and 
expansion are expected to be caused by peristaltic 
waves that propagate along the walls of ducts. In 
general during peristaltic action the fluid is pumped from 
lower pressure to higher pressure. Peristaltic transport 
occurs widely in the stomach, ureter, bile duct, small 
vessels etc. The principle of peristalsis is used by roller 
pumps for pumping fluids without being contaminated 
due to the contact with the pumping machinery. The 
initial work on peristalsis is done by using lab frame 
analysis. The important characteristics of peristaltic 
pumping namely trapping and reflux phenomenon are 
studied in detail by Shapiro et al [1] for the peristaltic 
flow of a viscous fluid through a tube and a channel. 

In physiological peristalsis the pumping fluid cannot 
always be treated as a Newtonian fluid. Kapur [2] 
suggested several Non- Newtonian models for 
physiological flows. He made theoretical investigations 
regarding blood as a casson and Hershel- Bulkley fluids. 
Blair and Spanner [3] reported that blood obeys casson 
model for moderate shear rate flows. Lew et al [4] 
reported that chyme is a Non Newtonian material having 
plastic like properties. In view of this Sreenadh et al [5] 



studied the effect of yield stress on peristaltic pumping 
of Non- Newtonian fluids in a channel. The Non 
Newtinian fluids are Bingham and Herschel- Bulkley 
fluids. Vajravelu et al [6, 7] made a detailed study on the 
effect of yield stress on peristaltic pumping of Herschel- 
Bulkley fluid in an inclined tube and a channel. All these 
investigations are confined to hydromagnetic study of a 
physiological fluid obeying some yield stress model. It is 
reported that some physiological fluids like blood are 
conducting fluids. Motivated by this MHD effects on 
peristaltic flow of a Bingham fluid in a channel with 
permeable walls is investigated in this chapter under 
long wavelength and low Reynolds number 
assumptions. The expressions for the velocity field in the 
plug flow and non-plug flow regions, the pressure rise in 
the channel and the volume flow rate are obtained. The 
effects of the magnetic field, yield stress and amplitude 
ratio on the pumping characteristics are discussed. 

II.MATHEMATICALFORMULATION 

Consider the MHD effects on peristaltic pumping of a 
Bingham fluid in a channel with permeable walls, under 
long wave length and low Reynolds number 
assumptions (Fig 1).The flow between the permeable 
walls is governed by Navier-Stokes equations whereas 
the flow in the permeable beds is according to Darcy's 
law. The channel is of half-width a. The region between 
y=0 and y=yp is called plug flow region. In the plug flow 
region | Tjr|<T[j in the region between y=ypand y=h, we 



have 



The wall deformation is given by 



H(X,t)=a+bsinH(X-ct) (1) 

Where b is the amplitude, A is the wave length and c is 
the wave speed. 

Under the assumptions that channel length is an integral 
multiple of the wave length A and the pressure difference 

across the ends of the tube is a constant, the flow 
becomes steady in the wave frame (x,y) moving with 
velocity c away from the fixed frame (X,Y). The 
transformation between these two frames is given by 
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X= X-ct; y= y; u(x,y)= U(X-ct,Y); 
Y); 



v(x,y)=V(X-ct, 



(-1+ it) + M C^Dy 



Where, e, = (^);, p =(-^ + m 



M 



M D 



Where (U, V) are velocity components in the laboratory 
frame and (u,v) are velocity components in the wave 
frame. In many physiological situations it is proved 
experimentally that the Reynolds number of the flow is 
very small. So, we assume that the wave length is 
infinite. So the flow is of Poiseuille type at each local 
cross-section. 

We introduce the following non- dimensional quantities 



We find the upper limit of plug flow region using the 
boundary condition that 



y/ =0 at >^ = >^ Q so we have yQ = — 

P (8) 



Also by using the condition z^^ = r,, at y=h 



h= a= ^ 



z 



We obtain P= 



h 



In the equations governing the motion (dropping the Hence — = — = r, 
bars) is ^ 



0<T<1 (9) 



("^o-y^J-^ '(?^, +1)=- 



dx 



0 = 



(2) 



(3) 



Where m = b 



(J 



The non- dimensional boundary conditions are 



0 



at y=0 



at y=0 



(4) 



(5) 



u 



= ¥=-l- 



^ D a du 



a 3 



at y=h 



(6) 



y 



Where y/ is the stream function, Da 
parameter and TqIs the yield stress. 



2 



, a= slip 



a 



equation (7) and using the relation (9) 

\\i in rvli in flrvxA/ rcin\r\m oo 



Taking y =yQ\x\ , 

we get the velocity in plug flow region as 



= —M C in M y M C 2C o s M y Q— ^ 



P 



(10) 



M 



Integrating the equations (7) and (10) and using the 
conditions y/^ =Oaty=0, y/-¥p ^y-y^- 

we get stream function as 



,—Py X 

y/ — C^osMy-\-C2^inMy + (, —)'> 



M 



2 



(11) 



The volume flux q through each cross- section on the 
wave frame is given by 



>'0 h 

q — ^ u ^dy -\- ^ u dy 



0 



yo 



P A n— M JiD 

(— -) + 



M 



2 



(12) 



M D 



SOLUTION OF THE PROBLEM 

Solving equation (2) subject to the boundary conditions 
(4) to (6) we obtain the velocity as 

p 

u —Xj/ =-M C ^SinM y + M C 2C osM y ^ (7) 

' M 



The instantaneous volume flow rate Q (x, t) in the 
laboratory frame between the centre line and the 
permeable wall is 

H h 

Q (X,t) = |t/ (X ,Y ,t)dY =|G + V}dy =q+h (13) 

0 0 

From equation (12), we have 
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d p 



3 

M D 



iq — s) 



2\q 



dx {k 2~ M hD 2) 



— M 



2 



(14) 



The average volume flow rate q over one wave period 
CT= 

^ of the peristaltic wave as 



Q =—\q{x ,t)dt =q +1 

T 



(15) 



Integrating the equation (14) with respect to x over one 
wave length, we get the pressure rise over one cycle of 
the wave as 



Kk^-M hD 2) 



Where h^x) = \+(l)sm27rx 

The time average flux at zero pressure rise (drop) is 
denoted by and the pressure rise required to 

produce zero average flow rate is denoted by A/?q so 
we have 



^P0 = -\ 



0 



(k^-M hD ^) 



X 



(17) 



The dimensionless friction force F at the wall across one 
wave length is given by 



F = 



dP ^ } fM 2(q —s) 

(k2-MhD2) 



\h(-—)dx =-\h{- 

^ dx *' 

Q ax Q 



-M ^}jx (18) 



I. DISCUSSION OF THE RESULT 



In order to see the effect of various pertinent parameters 

such as the yield stress C), amplitude ratio (^), Darcy 
number (Da), Hartmann number (M) and slip parameter 

(°) on pumping characteristics we have plotted Figures 2 
-11. 



The variation of pressure rise with time 
averaged flow rate is calculated from equation ((16)), for 
different amplitude ratios(e) and is shown in figure (2) 

for fixed cf= 2; Da=0.005; M=0.25; t=0.6. We observe 

that the higher the amplitude ratio, the greater the 
pressure rise against which the pump works. For a given 



, the flux Q for Bingham fluid depends on yield stress 
and it increases with increasing G5. 

From equation ((16)), we have calculated the 
pressure difference as a function of Q for different values 

of r is shown in figure (3) for fixed values a= 2; 
Da=0.005; M=0.25; e=0.6. It is observed that for a 
Bingham fluid, the peristaltic wave passing over the 
channel wall pumps against more pressure rise (ip) 

compared to Newtonian fluid. This type of behavior may 
be due to the presence of plug flow in Bingham fluid. 
Further, we observe that there is no difference in flux for 
Bingham fluid and Newtonian fluids for free pumping 
case (iip = 0). For a given iip the flux Q for a Bingham 

fluid depends on yield stress and it increases with 
increasing yield stress. From equation ((16)), we 

have calculated the pressure difference as a function of 
Q for different values of Darcy number Da, for fixed k = 

2; M=0.25; t=0.6; 0=0.6 and is shown in figure (4), We 

observe that for a given ip, the flux Q depends on Darcy 

number and it decreases with increasing Darcy number 
Da. For free pumping the flux Q is constant and it is 

independent of Da. 

The variation of pressure rise with time 
averaged flow rate is calculated from ((16)), for different 
values of « and is shown in figure (5) for fixed Da=0.005; 

M=0.25; T=0.6; 0=0.6. We observe that the longer the 

slip parameter, the greater the pressure rise against 
which the pump works. For a givenip, the flux Q 

increases with increasinga. For a given fluxij, the 

pressure difference increases with increasinga. 

The variation of pressure rise with time 
averaged flow rate is calculated from ((16)), for different 
values of M and is shown in figure (6) for fixed g= 2; 

Da=0.005; r=0.6; 0=0.6. We observe that the larger the 

Hartmann parameter (M), the smaller the pressure rise 
against which the pump works. For a givenip, the flux Q 

decreases with increasingiif. For a given flux^, the 

pressure difference decreases with increasingif . 

From equation (18), we have calculated the 
frictional force as a function of Q for different yield stress 

(t), amplitude ratio (S), Darcy number (Da), Hartmann 

number (M) and slip parameter (a) and it is observed 

that the frictional force F has the opposite behavior 
compared to pressure rise (Ap) and is depicted in figures 
(7-11). 
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Fig. 1 : Physical model for Bingham fluid in a channel with 

permeable walls. 
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Fig: 4.The variation of iij? \N\\hQ for the different values of 
Da for fixed e= 2; M=0.25; r=0.6; 0=0.6; 



FigT2.The variation of ij? with^ for the different values of 
0 for fixed cr = 2; Da=0.005; M=0.25; 1=0.6; 
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Fig: 5. The variation of Ap withfj for the different 
values of a for fixed Da=0.005; M=0.25; t=0.6; 0=0.6; 



: S.The variation of Ap with (J for the different^/alues^^ 
of T for fixed i£= 2; Da=0.005; M=0.25; p=0.6; 
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Fig: 6. The variation of A^^ with Q for the different values 
of M for fixed a = 2; Da=0.005; t=0.6; 0=0.6; 
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Fig: T.The variation of f^withKj for the different values of Fig: 10. The variation of f \N\\.hq for the different values 
0 for fixed e = 2; Da=0.005; M = 0.25; r=0.6; of cr. for fixed Da=0.005; M=0.25; 1=0.6; 0=0.6; 
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Fig: 8. The variation of \N\\hQ for the different values 
of X for fixed ^=2; Da=0.005; M=0.25; 0=0.6; 



Fig: 11. The variation of f with (J for the different values 
of jif for fixed ff= 2; Da=0.005; 1=0.6; 0=0.6; 

REFERENCES 

[1] Shapiro, A.H., Jaffrin, M.Y and Weinberg, S.L Peristaltic pumping 
with long wavelengths at low Reynolds number, J. Fluid IVIech. 
37(1969), 799-825. 

[2] Kapur, J.N. IVIathematical IVIodels in Biology and IVIedicine, Affiliated 
East-West Press Private Limited, New Delhi, 1985. 

[3] G.W.S. Blair, D.C. Spanner, An Introduction to Bioreheology, 
Elsevier, Amsterdam, 1974. 

[4] Lew, H.S., Fung, Y.C. and Lowenstein, C.B. Peristaltic carrying and 
mixing of Chyme in the small intestine ( An analysis of a mathematical 
model of peristaltic of the small intestine), J. Biomech., 4 (1971), 297- 
315. 

[5] Sreenadh, S., Narahari, IVI. and Ramesh Babu, V. Effect of yield 
stress on peristaltic pumping of Non-Newtonian fluids in a channel. 
International Symposium on advances in fluid IVIechanics, UGC-Centre 
for Advanced Studies in Fluid mechanics. Bangalore University, June- 
2004, 234-247. 

Fig :9. The variation of with (? for the different values of [6] Vajravelu, K., Sreenadh, S. and Ramesh Babu, v. Peristaltic 

transport of a Hershel-Bull^ley fluid in an inclined tube. Int. J. Non 
Da for fixed CJ = 2; M=0.25; 1=0.6; 0=0.6; Unear Mech. 40(2005) 83-90. 

[7] Vajravelu, K. Sreenadh, S. and Ramesh Babu, V. Peristaltic 
pumping of a Hershel-Bulkley fluid in a channel, Appl. IVIath. and 
Computation, 169(2005), 726-735. 




39 



www.ijitce.co.uk 



